Semi-classical model for the dephasing of a two-electron spin qubit coupled to a 

coherently evolving nuclear spin bath 
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We study electron spin decoherence in a two-electron double quantum dot due to the hyperfine 
interaction, under spin-echo conditions as studied in recent experiments. We develop a semi-classical 
model for the interaction between the electron and nuclear spins, in which the time-dependent 
Overhauser fields induced by the nuclear spins are treated as classical vector variables. Comparison 
of the model with experimentally-obtained echo signals allows us to quantify the contributions of 
various processes such as coherent Larmor precession and spin diffusion to the nuclear spin evolution. 



I. INTRODUCTION 

In recent years, electron spin qubits in solid-state quan- 
tum dots have emerged as promising candidates for the 
implementation of quantum information processing^"— . 
The confined electrons in these devices can be precisely 
manipulated using microwave frequency electromagnetic 
fields and/or nanosecond-scale pulses of nearby electro- 
static gates, while maintaining spin coherence over much 
longer times. The main source of decoherence in such 
qubits is the hyperfine interaction between the elec- 
tron spins and the nuclear spins of the host lattice^"—. 
Through this interaction, the nuclear spin bath produces 
a fluctuating effective Zeeman field on the electron spins. 
However, the timescale for evolution of this so-called 
"Overhauser field" is typically much longer than that 
required for manipulation of a single qubit. Therefore, 
dynamical decoupling techniquea^^'i^ based on fast con- 
trol of the qubit can be employed to partially eliminate 
decoherence due to the interaction with the nuclear spins. 
Recent experiments confirm that such techniques can be 
used to extend qubit coherence times by a few orders of 
magnitude, up to approximately 200 /xsiii^. 

Usually, interactions between a single electron spin 
and many weakly interacting spins (as described by the 
"central spin problem"—) can lead to complicated evo- 
lution of the joint quantum systemi^Sr— . Moreover, the 
nuclear spin bath may maintain coherence over a long 
time, which may in general result in coherent back-action 
on the electron spin. Indeed, using the Keldysh tech- 
nique and a non-trivial re-summation of diagrams, the 
authors of Rcfs. [32il33l analyzed this problem, and pre- 
dicted periodic collapses and revivals of the electron spin 
coherence over a specific range of external magnetic field 
strengths. This phenomenon was subsequently observed 
in two-electron spin echo measurements^^. 

Despite the apparent complexity of the system, there 
are several reasons why one might expect to find a more 
intuitive semi-classical description of the electron spin 
dynamics. First, due to the small nuclear Zeeman energy, 
the initial state of the nuclear spin bath is well described 
by an infinite temperature (completely random) state. 
Second, the state of the nuclear spin bath is not measured 
in the experiment, and the experimental outcome is an 
average over many runs. We note that the semi-classical 



approximation of the nuclear spin system has been used 
to describe a variety of other interesting phenomena in 
quantum dots^ii^iSai^. 

The aim of this paper is to demonstrate that a wide 
array of complex dynamical phenomena in two-electron 
spin echo measurements, such as those of Ref. [H, can be 
understood within the context of a semi-classical treat- 
ment of the nuclear spin bath. We first show that, within 
a simple model which neglects the effects of the Knight 
shift and the dipolar interaction between nuclear spins, 
the semi-classical treatment reproduces the expressions 
for the spin echo signal obtained in Refs. [s^lssl where a 
summation of diagrams in a perturbative quantum me- 
chanical treatment was used. We then present a more de- 
tailed microscopic model, which incorporates the Knight 
shift and the nuclear dipole-dipole interactions, as well 
as inhomogeneous hyperfine and Zeeman couplings. The 
semi-classical treatment for this model was sketched in 
the supplementary material of Ref. [isl Here we provide 
a systematic discussion of this semi-classical treatment, 
which relies on a low-order expansion of the qubit evolu- 
tion in the inverse of the number of nuclear spins and in 
the hyperfine coupling. Using this approach, we identify 
the relevant physical processes which govern nuclear spin 
evolution and electron spin coherence, as measured by 
the spin-echo signal. The excellent agreement between 
the model and the experimental data that was shown in 
Ref. [l^ lends additional justification to the approxima- 
tions underlying the semi-classical approach. 

The paper is organized as follows. In Sec. [H] we de- 
scribe the two-electron-spin qubit system and outline the 
main physical processes which govern the qubit dynam- 
ics. In Sec imi we review the Hahn spin-echo sequence. In 
Sec. II VI we describe the Hamiltonian of the two-electron- 
spin system and the nuclear spin bath. We then derive 
a simpler effective Hamiltonian through a perturbative 
treatment of the Overhauser field. In Sec. |V] we show 
that a semi-classical model based on this Hamiltonian re- 
produces the collapse and revival phenomenon. A more 
complete semi-classical approach is derived and justified 
in Sec. I VII There we include the effect of the dipolar 
coupling between nuclear spins and the back-action of 
the Knight field on the nuclear spin precession. Conclu- 
sions and discussion are presented in Sec. IVIII 
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II. THE TWO-ELECTRON-SPIN QUBIT 

We consider a qubit consisting of two electrons in a 
double quantum dot, in the regime where the two elec- 
trons arc well separated with one electron occupying each 
dot. A uniform in-plane magnetic field Bcxt = ^cxtz 
is applied along the z-axis, with i?oxt > 0. The qubit 
Hilbcrt space is spanned by the singlet spin state, \S) — 
(Iti) ^ lit)) I E^nd the triplet state with zero net spin 
projection on the z-axis, \To) = ^ (|tt) + lit))- In this 
notation, the two arrows represent the orientations of 
the electron spins in each of the two dots, measured rel- 
ative to the external field direction, z. Due to the large 
Zeeman splitting, this qubit subspace is energetically iso- 
lated from the two other two-electron triplet spin states 
|T+) = Itt) and |T„) = |||), in which both electron spins 
point parallel or anti-parallel to the direction of the mag- 
netic field. 

The energy splitting between the two qubit states \S) 
and \Tq) can be controlled on a fast timescale by rapidly 
tuning nearby electrostatic gates which modify the con- 
fining potential for the quantum dot, and hence control 
the shape of the two-electron wave function. For tunings 
where the two electrons are held far apart in the ground 
state, i.e. where the electrons are separated into differ- 
ent dots, the states \S) and |To) are degenerate. However, 
when the potential is tuned to favor partial double oc- 
cupation of one dot, the difference in orbital symmetry 
between \S) and \To) leads to an exchange energy split- 
ting J between them. 

In materials such as the commonly employed III-V 
compounds, the confined electron spins interact with a 
background of nuclear spins in the host lattice. This in- 
teraction is produced by the hyperfine coupling Hup ~ 
12d n ^d,nln 'Srf. Hcrc thc indcx n labels all nuclear spins, 
described by the operators {In}, d = L,R labels the elec- 
tron spins in thc left and right dots, described by the op- 
erators {Sd}, and the coupling constants {^d,n} depend 
on the local electron spin density, as will be described 
in detail below. Defining the nuclear (Overhauser) field 
operator in dot d as gVBBnuc,<i = J2n ^d.nln, we write 
the "effective" electron spin Zeeman Hamiltonian as 

Hpj = 9*1^3 (Si • Btot,L + Sfl • Btot.fl) , (1) 

with 

Btot.rf = Bcxt + Bnuc,d- (2) 

Here g* ss —0.4 is the electron effective g-factor in GaAs 
and fjiB is thc Bohr magneton. Equations ([T]) and ^ 
describe the Zeeman coupling in a system of two isolated 
electrons, where each electron is subjected to an effective 
field which is the vector sum of a uniform static external 
magnetic field Bcxt, and a local, operator- valued, Over- 
hauser field Bnuc,!i- 

We can gain extremely useful intuition about electron 
spin dynamics in this system by treating the operator- 
valued Overhauser fields Bnuc,L and Bnuci? as classical 



(time-dependent) vector variables. In typical GaAs dots, 
the Overhauser field is produced by a large number of 
nuclear spins in each dot, Nd ~ 10^. When all the nu- 
clear spins are polarized, the resulting effective Over- 
hauser field has a magnitude |Bnuc,ti| ~ 5 T. However, 
under experimental conditions, where thermal fluctua- 
tions randomize the directions of all nuclear spins, the 
typical values of |Bnuc,d| are reduced by a factor ^Nd, 
and are of order 1 niT. For strong enough external fields, 
Soxt > |Bnuc,d|, the net fields Btot,L and Btot,_R are 
nearly parallel to Bpxt- Under these conditions, the two- 
dimensional qubit subspace is only slightly perturbed by 
the misalignment of local fields, and remains energeti- 
cally isolated from the other two-electron spin states. To 

I B I 

leading order in ^""^"^ , the effect of the nuclear fields 
is simply to induce a Zeeman splitting between \\.'\) and 
Iti), proportional to the difference in z-projections of the 
effective fields in the two dots, ABfj^^ = B^^^ ^ — B^^^ ^. 

If we define a Bloch sphere for the qubit, whose poles 
on the z axis are the states Hf) and |ti), then the field 
ABj^^c points along the z axis. The states \S) and \Tq) lie 
on the X axis of this Bloch sphere. The splitting induced 
by AiJj^m, leads to oscillations between the qubit states 
\S) and I To), with a frequency proportional to |Ai3^^(.|- 
Such oscillations are polluted, however, by two sources 
of randomness. First, due to the fact that the nuclear 
state is random for an equilibrium nuclear spin bath, the 
magnitude of the initial nuclear field lAS^^cl; and hence 
the initial frequency of oscillations, is unknown. Further- 
more, due to internal dynamics of the bath itself, the nu- 
clear fields Bnuc,d(*) evolve in time. The resulting "spec- 
tral diffusion" of the qubit oscillation frequency leads to 
dephasing of the qubit oscillations^. Using an electron- 
spin-echo pulse, as explained below in Sec. lIIIi dephasing 
due to the unknown mean value of lS.B^^^{t) over some 
interval can be reversed. However, dccoherencc due to 
fluctuations of Ai?,^^^ (0 o'^i a timescale comparable to or 
shorter than the period between echo pulses in general 
cannot be eliminated in this way. 

For experiments involving weaker external magnetic 
fields and long enough evolution times, it is necessary 

IB I 

to go beyond the leading order in ' ^""•''i . Here we 
find that the transverse components of the Overhauser 
field, i?nucd(*)' crucially affect the qubit evolution in 
two primary ways. First, these transverse field compo- 
nents contribute to qubit decoherence by causing leakage 
of the electron spin state into the "non-qubit" subspace 
spanned by the states ||^) and |tt)- Second, the mag- 
nitude of the transverse part of the Overhauser field in- 
troduces a correction to the frequency of the \S) — \Tq) 
oscillations described above. Note that while in general 
a spin-echo pulse cannot reverse dephasing due to time- 
dependent fluctuations in the local fields, a partial or 
full recovery is possible if these fields vary periodically in 
time. Such a periodic time dependence, produced by the 
relative Larmor precession of different nuclear species, is 
the origin of thc collapse and revival phenomenon. 
An added complication in this moderate field regime 
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arises from the faet that, when we treat the Overhauser 
fields more properly as quantum mechanical operators, 

BnncA*)^ ^Lc,d,(0 and -Bnuc,d(0 do not commute. Con- 
sequently, at this order, the semi-classical approach to 
electron spin dephasing requires more justification. As 
we discuss below, such an approach is valid when the 
number of nuclear spins is large, and when the interac- 
tion between the electronic and nuclear spins is weak. 



III. THE SPIN-ECHO SEQUENCE 

With the above-described picture in mind, below we 
focus on the Hahn echo experiment in GaAs double quan- 
tum dots (see e.g. Ref. llSh . where each electron spin 
interacts weakly with N ^ 10^ nuclear spins. In such 
experiments, the two-electron state in the double quan- 
tum dot is initialized to the singlet ground state at large 
potential detuning, where both electrons reside in the 
right dot, 1(0,2)5'). Here the numbers in parentheses in- 
dicate the electron occupation numbers in the left and 
right dots, respectively, and the letter 5" indicates the 
two-electron (singlet) spin state. By rapidly tuning the 
potentials of nearby electrostatic gates, one of the elec- 
trons is transferred to the left dot within a time scale of 
approximately 1 ns. After this operation, the two sepa- 
rated electron spins evolve freely for a time ^ under the 
influence of the net local fields Btot,ci, given in Eq. ([2|). 
The gate potentials are then rapidly tuned to bring the 
electrons closer together. Here a substantial exchange 
energy splitting between \S) and |To) is maintained for a 
time corresponding to a "7r-phase" duration, which effec- 
tively leads to the the swapping of the states Hf) ^ Iti)- 
Then the gate voltages are rapidly tuned to separate the 
electrons, and the system is allowed to evolve over an- 
other interval of length ^ . 

At the end of the cycle, a spin readout procedure is 
performed. The gates are rapidly tuned to a large posi- 
tive potential bias, where the singlet ground state takes 
the "(0,2)" orbital configuration, while the orbital part of 
the triplet state remains of the "(1,1)" type due to Pauli 
exclusion. The charge configuration is then measured via 
a nearby charge sensitive detector. Due to the correla- 
tion between the orbital and spin degrees of freedom, the 
final spin state of the two electrons can be inferred from 
this charge measurement. Ignoring any imperfections of 
the measurement itself, we assume that the average of 
the charge detector signal taken over many runs depends 
linearly on the singlet return probability. 

The detector signal is averaged over a timcscalc which 
is long compared with all correlation times of the nuclear 
spin bath. Therefore we equate the averaged singlet re- 
turn probability for evolution duration r, Ps{t), with the 
average of single-run singlet return probabilities, taken 
over the equilibrium ensemble of initial nuclear spin con- 
figurations. Note that because the qubit is initialized in 
the singlet state, Ps{t) approaches 1 for very short evo- 
lution times T. For very long times, when coherence is 



lost and the qubit tends to an equal-probability classi- 
cal mixture of the states \S) and \Tq), Ps{t) approaches 
1/2. Therefore, by convention, we define the echo signal 
as 2Ps{t) — 1, which takes the value 1 for Ps(t) = 1 and 
for Ps{t) ~ 1/2. In this sense, the echo signal is used 
as a measure of electron spin coherence. 

In the recent experiment of Ref. [l^, the echo signal was 
observed to decay monotonically on a timescale of ap- 
proximately 30 IIS in high external magnetic fields above 
300 niT. At intermediate magnetic fields (120 - 300 niT) 
additional small fast oscillations were observed. At lower 
magnetic fields (50 — 120 mT) these oscillations evolved 
into a complete collapse of the echo signal at r sa 1 ^s, fol- 
lowed by a pattern of revivals and collapses on a timescale 
T sa 10 /iS. The collapse and revival pattern was attenu- 
ated by a decaying envelope over a timescale of r sa 30 fis. 

Our theoretical analysis of the echo experiments rests 
on a separation of timescales in the dynamics of the sys- 
tem. First, if one omits the 7r-pulse from the experimen- 
tal protocol, Ps{t) for the "free induction decay" decays 
to the value of 1/2 on the timescale T| 10 nsii^. This 
decay results from the uncertainty in the z component 
of the Overhauser field, AB^^^., which varies from run to 
run. In contrast, in a spin-echo measurement, the infiu- 
ence of a random, static Overhauser field AB^^^ on the 
final electron spin state is eliminated by the combination 
of the TT pulse and the two equal-length free evolution 
periods. If the Overhauser field were truly static, the 
electron spins would return to the state \S) at the end of 
the evolution. In a perfect measurement, for such a static 
Overhauser field, one would then obtain Ps{t) = 1, or 
an echo signal of value 1. Due to the time dependence 
of 'Bnuc,d{t), however, the echo signal typically decays to 
zero on a timescale of tens of microseconds. 

During the free evolution time while the electrons are 
separated, the system exhibits oscillations between \S) 
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and \To). Our crucial finding is that, for ' ^ 1^ the 

oscillations are well described in terms of tlie net accu- 
mulated phase determined by difference of magnitudes of 
the total effective fields on the two dots (we take h = 1): 

A$(t) = g>B f \ |Btot,L(t')l - |Btot.i?(i')l 1 dt'. (3) 



The magnitude of the total field is given by 
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tot.cZ 



|B^ . 



WP, (4) 



where |B^^j,^(<)| is the magnitude of the component 
of the Overhauser field in dot d which is perpendicu- 
lar to the external magnetic field. For -Bcxt not too 
large, |Btot,rf(^)l includes a significant contribution from 
l^nuc di^)\- Note that the time dependence of 'B^^^^ ^{t) is 
dominated by the relative Larmor precession of the three 
nuclear spins species, ^^Ga, ''^Ga and ''^As. Such rela- 
tive precession leads to a time-dependent modulation of 
|Btot,d| and causes a reduction of Ps{t) on a timescale 
of microseconds. In addition, random fluctuations of 
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-^nucdC*)' which arise due to interactions between nu- 
clear spins, lead to a reduction of Ps{t) on the same 
timescale. To account for all the above-mentioned ef- 
fects, which were observed in Ref. [l^ we thus include 
a combination of deterministic (Larmor precession) and 
stochastic processes in the evolution of the Overhauser 
fields in our semi-classical approach. 

IV. THE QUBIT HAMILTONIAN 

We begin our quantitative investigation by construct- 
ing the quantum Hamiltonian that describes the spin 
echo experiment. Although the experiment consists of 
several temporal stages, we focus on the free evolution 
periods, during which the Overhauser fields exert their 
main influence on the qubit evolution. We consider all 
other stages of the experiment, i.e. the singlet state ini- 
tialization at the first stage, the pi pulse, and the mea- 
surement of the final state, to be perfect. 

During the free evolution period while the electrons are 
well-separated, the qubit state evolves according to the 
effective Zeeman Hamiltonian 

Hpj = .9*M-B (Sl • Btot,L + ■ Btot.i?) , (5) 

presented above in Eq. ([Ij. The total effective field in 
dot d, Btot.d, is formed by the vector sum of the uniform 
external field Bcxt, and the Overhauser field 

g* fiB^nuc,d ~ ^^d,nln, ^d,n = ■A.a{n)\i'd,n\^ , (6) 
n 

where n is a label which indexes all of the nuclei. The 
parameter Aa{n) is the microscopic hyperfine coupling for 
nuclear spin species a(n), while the factor IV'd.nP weights 
the hyperfine coupling to nuclear spin n in dot d by the 
local electron density, and satisfies the normalization con- 
dition IV'd.nP = nc, where ric = 2 is the number of nu- 
clei per unit cell of the GaAs lattice. In regions where the 
electron density is substantial, the coupling to nuclei goes 
as |V'd,riP ^ where we define Nd = nl/J2n \i'd,n\'^ as 
the effective number of nuclei in dot d. For typical GaAs 
quantum dots, iV^ « (1 - 4) • 10*^. The index a = {1, 2, 3} 
runs over the three nuclear species ^^Ga, ^^Ga, and '^^As. 

To describe the evolution of the nuclear spin bath, we 
employ a Hamiltonian which includes Zeeman terms with 
a species- and site-dependent Larmor frequency for each 
nuclear spin, and the dipolar coupling between all pairs 
of nuclei: 

i/„uc = E^»^" + E^nx4^nM (7) 
n n.n' 

where n and n' run over all nuclei, and i and j label the 
Cartesian components of the nuclear spin operators. 

We wish to identify the main sources of decoherence 
for the two electron spin qubit which arise from the com- 
bined evolution under the Hamiltonian in Eqs. ([T]) and 



Naturally, we shall use 1/Nd as a small parameter. 
In addition, we note that the lowest external field used 
in the experiment in Ref. [3 (50 mT) was more than an 
order of magnitude bigger than the typical magnitude of 
the Overhauser field. Therefore we will proceed to study 
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decoherence effects as an expansion in ' ^""'^'i ^ and to 
leading order in 1/Nd- 

First, note that if one replaces the operators 'Bnvic,L(R) 
in Eq. © by classical vectors with magnitudes much 
smaller than Boxt, then, similar to the case when the 
Overhauser field was absent, the system described by 
Hamiltonian ([T]) possesses a two-dimensional subspace 
which is energetically well-separated from the remaining 
two levels. This new qubit subspace is spanned by the 
states Itfti) ® linn) and \in^) ® Wur), with eigenvalues 
±5MBff* (|Btot,L| - |Btot,fl|)- Here the up and down ar- 
rows indicate the projections of the electron spins onto 
the quantization axes fi^ and hn parallel to the total 
fields Btot,L(_R) in each dot. 

Deviations of the directions of and fiji from the 
z-axis arise from the Overhauser field components per- 
pendicular to the applied field, B;j-^^ ^. Due to evolution 
of the nuclear spin bath (primarily due to Larmor preces- 
sion of the nuclear spins), the fields B^^^^^^ (i(^)' hence 
fiL and fiR, are slowly modulated in time. However, be- 
cause the frequencies of such modulations are typically 
two orders of magnitude smaller than the value of the 
energy gap between the two instantaneous eigenstates, 
IS.E = ^b9* (|Btot,L| - |Btot,fl|), we assume that each 
electron spin adiabatically follows its local, slowly vary- 

^ IB I 

ing, quantization axis fiL or ur. For small , the 

main effect of the nuclear field is thus to modulate the 
magnitude of the total field |Btot.(i(i)|, Eq. (j4)), and hence 
to modify the dynamical phase / AE{t)dt accumulated 
between the two eigenstates over the free evolution pe- 
riod. We therefore ignore changes in the directions of the 
quantization axes in each dot, and describe the evolution 
of the system by using the effective Hamiltonian 

HeL. = g*f^B {SI |Btot,L| + |Btot,ii|) . (8) 

For given classical values of Btot.L and Btot,i?: this 
Hamiltonian preserves the instantaneous eigenvalues of 
the original Hamiltonian H^i, Eq. ([T]). Expanding |Btot,d| 
in Eq. (|4]) in the small parameter ^^"°^|'^ , the Hamilto- 
nian Hei z becomes 

i?eZ..«5>B E f^nuM + ^m^l^J- (9) 
d=L.R \ I cxtl y 

The term proportional to B^-^^^ ^ is the zcroth-order con- 

I B I 

tribution in and may give rise to a magnetic- 

field-independent contribution to the electron spin deco- 

IB^ 1^ 

herence. The term proportional to i'b^'^i is the first 

I B I 

order correction in |^"''^| , and is responsible for the 
interesting collapsc-and-rcvival behavior that we study 
below. 
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It should be noted that in writing Eqs. ([8]) and ([9]), 
we have ignored effects arising from the relative angle 
Q \^^<^,d\ ijg^^ggjj ^j-^g quantization axes in the two 
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dots, which enter at order 



|B„ 



First, the misalign- 



ment of axes may cause unwanted transitions to the \T±) 
states when initializing from the singlet state, or during 
free evolution under Hamiltonian i/gj. However, these 
effects lead to a reduction of Ps on the order of 9^ . Sec- 
ond, we neglect any possible geometric phases which may 
accompany the dynamical phase accumulated while the 
electron spins adiabatically follow the local quantization 
axes in their separate dots. Such phases are proportional 
to the areas of closed loops swept out by and hn dur- 
ing their evolution, and for small 9 are also proportional 

to 

Although we obtained Eq. ^ by treating the nuclear 
spin operators as classical vectors in a Taylor series ap- 
proximation to Eq. an identical expression for ifgZ.z 
was formally derived in Ref. [s^ from the full quantum 
Hamiltonian, Eq. ([Ij, using a SchriefFer- Wolff transfor- 
mation. The classical argument above thus provides a 
simple intuitive explanation for the formal perturbative 
derivation. Hereafter, unless otherwise specified, we treat 
the fields in Eq. (j9]) as quantum operators with appro- 
priate commutation relations. 

To complete the description of the problem, we now 
discuss the 7r-pulses employed in the Hahn echo sequence. 
These pulses are achieved by applying a time-dependent 
perturbation i?7r(i), which adds to the system's total 
Hamiltonian, H = H^lz + Hnuc + H-Kit)- We assume 
HT^{t) is only nonzero over narrow intervals which arc 
short compared to all timescales relevant for evolution 
under Hf-i^z and -ffnuc- Rather than specifying a detailed 
time-dependent protocol for (t) , we define H^^ implic- 
itly in terms of its effect on the electron spin operators 

CZ. 

f U!^HAt')dt'^ ^^-r [e-'^o = c(i)5J, (10) 

where T (T) is the (reversed) time-ordering operator. 
In writing Eq. pU)) . we assume that H.,j{t) acts only 
within the two-dimensional qubit subspace. We shall 
consider perfect 7r-pulses, for which the "echo function" 
c{t) switches between 1 and —1 over the short duration 
of the pulse. For simplicity we consider the pulses to be 
instantaneous, and for the Hahn echo sequence write 



c{t) = e (r/2 - t) - e (t - r/2) , 



(11) 



Using Eqs. pll| and pT|) we switch to an interaction 
picture with respect to HT^it), where 



Sl{t)^c{t)Sl 



(12) 



becomes 



d=LM 



nuc.fi I 
2 1 -Bcxt I 



c(t)5J. 



and where we employ the notation 5| = {t = 0). The 
interaction-picture time-dependent Hamiltonian H^i^zit) 



(13) 

Equation p^ . with Eq. ([7]), will serve as the starting 
point for our analysis of decoherence in the spin-echo 
experiment. 



V. SEMI-CLASSICAL MODEL FOR THE 
REVIVALS 

We now show that a simple semi-classical approach in 
which we treat the Overhauser field operators as clas- 
sical vectors can reproduce the electron spin coherence 
collapse and revival effect predicted in Ref. [s^ and ob- 
served in Ref. [isl In the next section, we will provide a 
systematic derivation and justification of this approach, 
starting from the full quantum description. 

In this section, we treat the Overhauser field in each 
dot d in Eq. as a sum of three classical vectors, 

Bnuc,d(i) = X]q=i BQ.d(i), where a indexes the three nu- 
clear spin species. We assume that the magnitudes of 
the species-dependent fields {|Bq,^|}, and their z compo- 
nents {B^j}^ are random but constant throughout the 
evolution. The time-dependence of Bnuc,d(i) within each 
run arises solely from the Larmor precession of the trans- 
verse nuclear spin components. Explicitly, we neglect the 
nuclear dipole-dipole interaction, and the influence of the 
Knight shift on the nuclear Larmor precession. We as- 
sume that all nuclei of the same species precess at a single 
Larmor angular velocity, uja = 7a-Bext- 

The echo signal Ps{t) is obtained by averaging over 
many experimental runs. Thus we must average the re- 
sults of electron spin dynamics against the distribution 
of initial states of the nuclear spins. Due to the large 
number of nuclear spins, N sa 10^, the initial values of 
the components of each vector 'Bandit = 0) are Gaussian 
distributed with zero mean and a standard deviation ba^d 
of order 1 mT (sec calculations below). 

The model in Eq. (|13p . under the assumptions above, 
is sufficient to produce the collapse and revival ef- 
fect in Ps{t), and further provides an intuitive semi- 
classical picture in which to understand the phenomenon. 
However, because we neglect the time dependence of 
-^nucrf(*)' ^^'^ effects of the Knight shift and other 
dcphasing mechanisms of the nuclear Larmor precession, 
this model does not capture the decaying envelope ob- 
served in the experiment of Ref. [l^. These issues will be 
addressed in detail in Scc lVIl 

We now calculate Ps{t), using the singlet initial state 
m = 0)) = \S) = ^ (I n) - I it))- For any given set of 
initial values of the (18 total) components of the six clas- 
sical vectors {'Bandit = 0)}, the Hamiltonian in Eq. 
generates a pure quantum evolution which after an evo- 
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lution time t = t yields 

p-iA*(r)/2 

mr)) = 



V2 



,iA<I.(T) 



(14) 



The relative phase A$(t) = $i — is related to the 
difference of the magnitudes of the total effective fields, 
|Btot,d(i)|, hi the two dots. Within the approximation of 
|Btot, ci(i) I used to write Eq. p^ . we obtain: 



$d(T) 



g MB 

2\B, 



cxt I 



IBic.dWI'cWdt. (15) 



Note that the phase $c; is determined solely by the dy- 
namics of a single isolated electron in dot d. Below we 
will use this fact to relate the decoherence of the two- 
electron singlet-triplet qubit to that of a single electron 
spin in a quantum dot. 

For the final state \iP{t)) in Eq. (|14p . the singlet return 
probability is given by 



\{S\^{t))\' 



1 



cos(A<I>). 



(16) 



The ensemble-averaged singlet return probability Ps{t) 
is found by averaging the result for a single run, Eq. (|16p , 
with respect to the distribution of initial magnitudes and 
directions of the six vectors {'Ba.d{t = 0)}. Note that be- 
cause cos(A$) — Re [e**^e~'*^] , and because the Over- 
hauser field configurations in the two dots are indepen- 
dent, we can average over e'*^ and e**^ independently. 



1 1 



){e- 



(17) 



To perform the averaging, we calculate ^^(t) us- 
ing Eq. ([T5|) . with c(t) given by Eq. ((TT|) . and with 
the classical evolution for Bj|^j,j^(t) resulting from the 
free precession of the underlying nuclear spins. We 
write \B^^,Jt)\^ = Ec.fsb^^t^Jt) in terms of the 

complex- valued fields ba^d = d~^^^^a d- ^^^^ decom- 
position, the time evolution of the Overhauser field due 
to Larmor precession is given by ba,d{t) = &Q,d(0)e"^'»*. 
Each phase ^d is then given by 
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215, 



cxt I Jo 



c{t)Y,b^^d{t)hl^d{t)dt 



(18) 



Kd{0)bl^d{^)— sin^ {uJc.pT/ A) 



a, 13 



where uja/s — uJa — ujp. 

Next we must integrate over all initial conditions, i.e. 
over the initial magnitudes and phases of the three fields 
{&c(,d(0)} in each dot d. For this purpose we express 
the initial conditions as 6^,^(0) = ba^dZa, where each 
Za = Xa + iya is a dimensionless complex variable. The 
quantity fo^d is the root-mean-squarcd (rms) value of 
each component of the transverse field for species a in 
dot d, 



(19) 



|(/q + 1)/q, = |. Here Ua is the average 



with Ga 

number of nuclei of species a, per unit cell (see Ref. |3 
This substitution gives 



2_^Tap.d- 

a,0 



(20) 



with 

Tal3,d 



IB, 



ILOap 



sin" (Wa/3T/4) . (21) 



We now carry out an ensemble average over the 
initial conditions by treating all components of the 
{zq,} as independent Gaussian-distributed random vari- 
ables with zero mean and unit variance, i.e. with 
probability density function p({zq.'}) J^^^j^ c?a;Q,(iya = 



■nLiexp(- 



)dxadya: 



a 



(22) 



Here the parameters {Xa,d} are the eigenvalues of the 
T-matrix for dot d. 

The 3x3 Hermitian matrix T in Eq. ([21]) corre- 
sponds to that of Rcf. 32. Because T is Hermitian 
and is similar to an antisymmetric matrix, it has one 
zero eigenvalue, Ai^^ = 0, and a pair of eigenvalues 

Inserting these eigenvalues 



A 



2(3), d - ±\/Sq>^ \-LaP,d\ 

into Eq. ([2^ . and using Eq. (|2ip for T^^.d, we obtain: 



E 

a>l3 



Ag* tJ-Bba4bfj^d\ sin'' (i:Ja^r/4) 



UJ 



(23) 

Thus we see that the semi-classical model used in this sec- 
tion reproduces the result of Ref. [sH for the decay of the 
spin echo signal in a single quantum dot. The echo sig- 
nal shows oscillations with amplitude (^-^]^~p;3~f 
which develop into the complete collapses and revivals 
at low magnetic fields. Note that the expression for 
the phase in Eq. (|18p includes terms which are bilin- 
ear combinations of Gaussian variables, and so are not 
Gaussian themselves. Therefore the decay in the interval 
T ^ I/wq/j behaves like an inverse polynomial, rather 
than the form e~'^°"''*''^* expected for spectral diffusion 
(see below and Refs. [35l - [37t) . 

We now return to computing the echo signal in the 
double dot system, Eq. ([TT]). First, note that (e"**'') in 
Eq. is strictly real. This fact is a consequence of 
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antisymmetry of the echo function around the time r/2. 
Thus we can drop the "Re" from Eq. ([T7)) and write the 
echo signal as 



2Ps 



1 



(24) 



expansion i^^^^^^i 
which involve nuclear 



The expressions in Eqs. ([22)) and ([23)) . which describe 
the dephasing of a single electron spin in an isolated 
quantum dot, were derived previously from a fully quan- 
tum mechanical treatment in Ref. l32l A key element 
of the derivation in that work was the vanishing of the 
contribution of the commutator [I'^ between nuclear 
spin operators in the low order perturbation expansion 
of the evolution operator. The vanishing commutator is 
indicative of classical behavior, and further motivates our 
classical treatment of the nuclear evolution. 

To conclude this section, we show that the clas- 
sical treatment provides an intuitive explanation for 
the collapses and revivals. The total effective elec- 
tron Zeeman field in each dot d, |Boxt + Bnuc.dl, de- 
pends on the square of the transverse Overhauscr field, 
l-^nucdP- '^^'^ terms in the expansion IB-'- 
V ( + B^ B^ 

spins of two different species, a ^ a\ oscillate at the 
relative Larmor angular velocity uja ~ '^a' ■ As a result, 
the magnitude of the total field, IBtot.dl, and hence the 
splitting between electron spin energy levels, oscillates 
as a function of time. These oscillations determine the 
phase accumulation during each run of the experiment. 
In this simple semi-classical treatment the nuclear spin 
evolution is not affected by the electron spin, so that the 
phase between electron spin components remains well- 
defined during each run. However, averaging over the 
ensemble of nuclear spin initial states amounts to averag- 
ing over the phases and amplitudes of these oscillations, 
and causes the collapse of the echo signal. However, if 
the free evolution time t/2 is simultaneously an integer 
multiple of each of the three relative Larmor periods, 
then the contribution of those oscillations vanishes inde- 
pendently of the initial nuclear state, and a revival of 
electron spin coherence is observed. Note that it is a for- 
tunate coincidence in GaAs that the Larmor frequencies 
of the three species are nearly equidistant, which, to a 
good approximation, allows the commensurability condi- 
tion to be easily fulfilled simultaneously for all three pairs 
of nuclear species. 



Eq. ([13)) . We also now reintroduce the nuclear dipolar in- 
teraction, see Eq. ([7]); and allow for local variations of the 
nuclear Larmor frequencies, which were omitted from the 
simple model in section|V| At each step in the derivation, 
we will justify the approximations needed to arrive at the 
semi-classical model with arguments from first principles, 
or with the help of experimental data. Along the way, 
we will also incorporate several important features of the 
dynamics, such as the time dependence of i3^„(,(f), back- 
action of the electron spin evolution on the nuclear state, 
and nuclear spin dephasing, which were ignored in the 
heuristic treatment in the previous section. 



A. Quantum expression for the echo signal 

Our aim is to derive an expression for the echo signal at 
the end of the full Hahn echo sequence, via the calculation 
of the singlet return probability 



Ps{t) = ^Tr„uc 



{S\U\t)\S){S\U{t)\S)], (25) 



where U{t) = T-g-* JJ" (^ei,=(t)+^^nuc)dt jg ^^^^ interaction 
picture evolution operator of the joint electron-nuclear 
spin system, with T representing the time-ordering oper- 
ator. Here we have assumed that the electron spins are 
initialized to the singlet state \S). The trace in Eq. ([25)) 
is taken over all nuclear spin states, with Z^o = Trnuc[l] 
representing the partition function for an infinite tem- 
perature (completely random) nuclear spin state. 

To separate the time evolution due to the static part 
of -B^uc (which dominates the evolution for t < r, but 
has no influence onU{t = t) due to the spin-echo), we 
introduce the zero-order Hamiltonian 

Ho{t)=gyB E B,U<i'^W^d+E^fe4^ (26) 

d=L.R k 

and a corresponding zero-order evolution operator 
f7o(t) — e~^^o ^o(t )dt _ Note that due to antisymmetry 
of the echo function, /J" c{t')dt' = 0, the evolution oper- 
ator Uo{t = r) at the end of the full sequence does not 
depend on the electron spin operators and simply rotates 
all nuclear spins about the z-axis. 

The full evolution operator U (r) can be rewritten as 

U{t) = Uo{t) ■ re"''''o'''**^=i[-f^^-''(*)^^'=(*)+-f^°-''(*)l, (27) 



VI. DERIVATION OF THE SEMI-CLASSICAL 
APPROACH 



In this section we present a systematic derivation of the 
semi-classical treatment of electron spin dynamics pre- 
sented above. Keeping in mind the discussion surround- 
ing Eq. ([9]), we now restore the quantum nature of the 
electron spin operators and begin with the interaction- 
picture quantum Hamiltonian H = Hei^z{t) + H^uc, 



where the time-dependent operators in the exponent 



are given by Hj_,d{t) = U^it) 



•9 MB 2|B„,| 



^n.n'^d n,n' n n' 



Uo{t) and 
Uo{t), see Eq. 1^. 



These terms describe the evolution of the Overhauser 
fields due to Larmor precession and the dipole-dipole in- 
teraction between nuclear spins in dot d, respectively. We 
assume that the two dots are well separated, such that 
the inter-dot dipolar coupling can be neglected. This 
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approximation is not essential for the derivation ahead, 
however, and the existence of a smah inter-dot couphng 
would not significantly affect the final result. 

In order to evaluate Eq. (PS)) . we decompose the evo- 
lution operator U (r) into four separate pieces. Because 
operators describing spins in different dots commute, the 
exponentials in Eq. ([27| can be factored by dot index 
d. Therefore the evolution operator can be written as 
U{t) ~ Ul{t)Ur{t), where Ui^(^ji^{t) only involves spin 
operators in the left (right) dot. In addition, because 
the z-projection of electron spin in each dot is conserved 
by C/(t), see Eq. ^7}, we can separate the evolution by 
introducing projectors onto the product states | f^) and 
lit): 

U{t) = UL+{T)UR-{r) in> {n\+UL-{r)UR+{r) lit) Utl , 

(28) 

where Uda is a unitary operator acting only on the nu- 
clear spins in dot d, with the electron spin taken to be in 
the eigenstate | t) for tr = +, or | J,) for ct = — . In- 
serting Eq. psp into Eq. <\2b\ . and performing the trace 
over nuclear spin states separately for the two dots, we 
find that the echo signal is given by [c.f. Eq. (p4|) ] 



2Fs - 1 = Re {Cl ■ Cr) , 



(29) 



where the superscripts -I- and — indicate projection onto 
the I t) and | l) electron spin states. The projection is 
needed because H±{t) and Ho{t) depend on Uo{t), which 
ior t ^ T depends on S^. 

Intuitively, the two terms involving H±{t) and Ho{t) 
lead to suppression of the echo signal, each through a 
different physical mechanism. In the spirit of the semi- 
classical approximation in Sec. [Vj where the Overhauser 
field operators were treated as classical vector variables, 
H±{t) causes dephasing through the time-dependence of 
l-^nuc(*)l generated by nuclear Larmor precession, while 
Hoit) causes decoherence through fluctuations of B^^^{t) 
generated by dipolar-interaction-mediated spin diffusion. 

The semi-classical approach relies on three main ap- 
proximations, which involve neglecting various commu- 
tators of the form [B^„^(t), BJ„^(i')]. When the num- 
ber of nuclei Nd is large, these commutators scale as 
\^-mic{i)\l Nd- Below we provide a physical motivation 
for each such approximation, and discuss the associated 
range of validity. 

First, we neglect the commutator between H±(t) 
and Ho{t'), and factor the exponentials appearing in 
Eq. dSJ): 



where 
Cd 



-Tr„ 



C/]^(T)[/rf_(T)l , d = L,i? (30) 



The echo signal, Eq. ([29l) . is a product of two similar 
averages which arc taken over the different sets of nuclear 
spins in the two dots. In the remainder of this section we 
focus on the behavior of the average within a single dot. 



— Tr 

2 -LInuc 



t/i(r)C/_(r) 



This 



which we denote as C 

allows us to simplify all formulas by suppressing the dot 
index d. At the end of subsection IVI CI we will return 
to the two-electron double-dot echo signal, including the 
combined effects of dephasing in each of the two dots. 



B. The semi-classical approximation of separating 
the dynamics of the spin diffusion and transverse 
Overhauser field 

In this subsection we introduce a semi-classical treat- 
ment in which the evolution operators in Eq. ([29| are fac- 
torized into contributions depending separately on H± (t) 
and Hd {t) . We start by using Eq. ([27l) and the fact that 
Uq{t) involves only nuclear spin operators to write the 
evolution operator in C as 



uI{t)U^{t) = T 



^^gdt[^H+{t)c{t}+H+{t)] 



(31) 



J- \^-\f<; dt'[-^Hl{t')c{t')+H^{t')] 



with 



(32) 



(33) 



The commutator [H±{t), HD{t')] which was set to zero 
in writing Eq. ([32]) leads to nuclear spin dephasing, i.e. 
decay of the correlator {Il{t)ll{t')) , due to dipolar flip- 
flop events. In our model, we account for this effect 
phenomenologically by introducing random site-to-site 
small variations of the nuclear Larmor frequencies (see 
Sec. IVI C]) . In doing so, we assume that intrinsic and ex- 
trinsic dephasing of nuclear spins affect the electron spin 
dynamics in the same way. 

In addition, the commutators neglected above may 
lead to enhanced nuclear spin diffusion through the com- 
bination of dipole-dipole and hyperfine-mediated spin 
flips. The validity of this approximation thus depends 
on material and system parameters. However, it can be 
checked by comparison with experimental data, as dis- 
cussed in more detail in section. I VIII 

The second approximation is to split the average of 
the product of operators involving Uoa and Uj_a into a 
product of averages involving Uoa and Uj_a separately: 



9 



-Tr„ 



-Tr„ 



-Tr, 



(34) 



r 



The correlations which are neglected by splitting this av- 
erage are related to correlations of the longitudinal and 
transverse components of the Overhauser fields, B^^(,(t) 
and B^jj(.(i). These correlations are contributed by ex- 
pectation values of at least four operators of the same 
nuclear spin n, e.g. (/^) (/^) . On the other hand, the 
leading contribution to the split-average comes from the 
expectation values of only two operators of the same nu- 
clear spin, e.g. (/^)^, and therefore is roughly Nd times 
bigger, due to the fact that it contains a sum over at least 
Nd times as many terms. 



C. The semi-classical approximation of averaging 
over all nuclear spin states 

Through Eq. ((34| . the dephasing in a single dot can 
be approximately expressed as a product of two sepa- 
rate averages, one involving Udo- and the other involv- 
ing U±a-- The average over Uoa- can be associated with 
decoherence due to spectral diffusion caused by dipole- 
dipole mediated nuclear spin diffusion. This process was 
analyzed in Refs. [ssI - IstI . and was shown to result in an 
echo decay factor cxp{~ (t/tsd)'^)- Using this result and 
Eq. (1211), we rewrite Eq. dSHl) as 



C(t) 



Trn 



C/l+(r)[/^„(r)l , (35) 



The spectral diffusion time tsd depends on the details of 
the quantum dots, and weakly depends on the magnetic 
field strength. According to Refs. [35ll36l tsd is of order 
10 microseconds, which is comparable to the duration of 
the experiment. Thus, the decay factor is likely to be im- 
portant for a detailed fit to the experiments, and will be 
dominant at high magnetic fields. However, as we show 
below, at low enough magnetic field the dephasing asso- 
ciated with C/_Lcr, i.e. with the evolution of the transverse 
components of the Overhauser fields, becomes dominant. 

We now turn to calculate the the dephasing associated 
with U±cr semi-classically. Making use of the large num- 
ber of nuclei involved, Nd > 10^, the trace in Eq. (PSI) can 
be cast into a form which closely resembles the classical 
average described in Sec. IVl We first combine together 
large groups of nuclei with similar couplings to form a 
set of "giant" collective spins. Then, by evaluating the 
trace in a basis of coherent states with well defined ori- 
entations of these giant spins, we find that their induced 
Overhauser fields effectively act as classical variables like 
those introduced "by hand" in Sec. IVl However, there 
are two important differences in this more refined treat- 
ment. First, whereas in Sec. IVl all nuclear spins within 



each species were forced to precess at a fixed Larmor fre- 
quency, we now account for the electron-spin-dependent 
shift (i.e. the Knight shift) of the nuclear precession rate 
which is inherent in Eqs. Q and (|26p . Second, we now 
account for inhomogcncity in the system, both in the 
Knight field terms and in the Larmor frequencies of all 
nuclear spins. Site-to-site variations of the Larmor fre- 
quencies Un, see Eq. ©, phenomenologically account for 
nuclear spin dephasing due to, e.g. local quadrupole mo- 
ments. We assume that the inhomogeneities are weak, 
such that the differences between the mean Larmor an- 
gular velocities of the three species, LOa ~ Ta^cxt, are 
much larger than the widths of the distributions for each 
one. 

To account for inhomogeneities in the system, we di- 
vide the nuclei in each dot into K groups labeled by 
the index k — l,...,K, with each group k containing 
nuclei of the same species a{k). These Nk are picked 
from the Nd nuclei in dot d such that all members of the 
group have nearly the same hyperfine coupling Ak and 
feel nearly the same phenomenological local shift in mag- 
netic field, SBk- The Larmor angular velocity for nuclear 



spins in group k is given by uj^ 



\ac(t)Ak where 



5ujk = "1a{k)^Bk- Note that the sign of the Knight shift 
±i^fe depends on the state of the electron spin in the 
dot at time t, and is therefore proportional to ac{t). 

To calculate C, the single dot contribution to the 
echo signal, we need to choose a basis for the nu- 
clear spin Hilbcrt space. For each group k, we 
form a collective spin from all of the nuclear spins in 
the group, and consider "giant spin" states of well- 
defined total angular momentum Ik and orientation 
ftfc = (sin6'feCOS(/9fe, sin6'fcsin(^fe, cos9k). |Ife|^|4,nfe) = 
Ikih + l)|/fc,nfc), (fife • Ife)|4,nfc) = /fc|4,nfc), where 0k 
and (fk are the polar and azimuthal angles of giant spin 
k, respectively. A wave function in the Hilbcrt space of 
all nuclear spins is written as j^*) = <S>k=i l^k^'^^k)- The 
trace is then performed by summing over all possible as- 
signments of the values {Ik} with appropriate weights^, 
and integrating over all possible directions {iik}- The 
length of each giant spin k can vary from to ^Nk, 
but the vast majority of giant spin states have lengths 

of order Ik ~ \/ Nk- The relative quantum uncertainty 
in the transverse components of the giant spin coher- 
ent states, Aiyil, scales as l/\/Tk = l/^/Wk. As we 
shell see below, for analyzing the behavior near the re- 
vivals peaks it is sufficient to divide the nuclei into about 
K ^ \Q groups, which leaves as many as Nk ~ 10^ nuclei 
in each group. 

In the coherent state basis, the semi-classical 
approximation converts the expectation values 
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{'i'\U]_j^^{T)U±d-{T) \'^) , which arise in the evalua- tion of the trace in Eg. (1351) . into "classical" expressions 



(36) 



where IB^^^J^ ^(t)|^; is determined by the expectation values I^^'"^\t) ~ ( J^, ftfc | 7^'^''''' (i) | /fe, fi^ ) of the giant spin 
components: 



1 



4:g*fj,BBc 



^ A,A„ (ii-it) ■ qnt) + 17 {t) ■ lint) 

k,k' 



The superscript a indicates that the time dependence of 
the operators is determined by the evolution with respect 
to Ho, Eq. which depends on the electron spin state 
a. 

The validity of the approximation, Eq. p6p . is dis- 
cussed in Sec. IVI El below. After making the approxima- 
tion, however, the trace reduces to an integration over all 

possible expectation values l7^^\t) of the components 
of the K giant spins. We perform the integrations us- 
ing a similar method to that used in Section |Vl First 
we write the classical field \B^^^{t)\li = Y.k,iK^t)^l* i^) 
as a sum of products of complex variables, where 6^ = 

I 



^fe {iTi^) + ilk" it))- Each bl{t) evolves according to 

bl{t) = l)f^(^Oy['^.ik)t+Su:kt+aAkf^cit')dt']^ ygjj^g ^-j^g 

proximation in Eq. (|36|) to evaluate the trace in Eq. psp , 
we have 



C(T)«C,e(r) = e-(^/^--''(e'*) 



(37) 



where the angled brackets indicate a "classical" averaging 
over the complex random variables &^(0), and the phase 
$ is given by the integrals in the exponents in Eq. 



<i>(r) = 



9 fJ-B 



cr=±l 



^ / dtY,but)br{t) 



k.l 



■Y,bk{0)bi{0) f dtc{t)exp 
kJ a=±i-'o 



i{ujki + 5ujki)t + icrAki / dt'cit') 
Jo 



(38) 



with LOkl = U!a{k) - ^^a{l), SuJkl 

Ak - Ai. 



Sujk 



Sllii, and Am 



Assuming that Nk is large for every group fc, we per- 
form the average over all initial coherent nuclear spin 
states by writing bk{0) — bkZk in terms of a collection of 
independent, complex Gaussian random variables {zfc} 

with unit variance. Here g* pLgbk - 



A^/caa(fe)/2Afc is 
the rms value of each component of the Overhauser field 
associated with group k in the dot, similar to the param- 
eters {6a} appearing in Eq. pip . In this representation, 
Eq. ((55)) becomes 



ZkZi 



(39) 



k,l 



with 



Tm (r) 



X 4 



ig* fiBbihjuJki + SuJki) 

2|i?ext| 

C0s(Amt/2) - cos [{uJkl + SLL!kl)T/2] 



(40) 



- All 



Performing the Gaussian average over {z^}, we obtain 



k.l 



= j^_dxk.dyk.p{\zk.?)e^v\-^Y.Tkl^ 
1 



+ «A„i (r) ' 



(41) 



where the parameters {A^} are the eigenvalues of the 
M X M Hermitian matrix with elements {T^i}- 
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Finally, we calculate the spin-echo signal of the two- 
electron double quantum dot singlet-triplet qubit, incor- 
porating the dcphasing due to both two dots, Eq. ([29| . 
using Eqs. and (gT]) for each dot L,R: 



and I comprised of the same nuclear species), we find 



2P, - 1 



n (niT^ 



(42) 



d=L,R \ m 



Here tsd is the effective spectral diffusion timescale for 
the two electron system, f^^ = t^^ ^ + t^^ ^. For non- 
identical dots, the eigenvalues Xm,d are generally different 
due to the differing number of nuclei Nd, and their associ- 
ated distributions of coupling constants {Ak} and {w^}. 
All of these parameters affect the grouping of spins into 
"giant spins," and the matrix elements Tik-, calculated 
according to Eq. 

Equation (|42|) is the main result of the paper. This re- 
sult shows that, within the semi-classical approach, the 
echo signal 2Ps — 1 in the Hahn echo experiment can 
be understood in terms of a decay envelope arising from 
spectral diffusion, along with an additional factor which 
arises from the relative precession of different nuclear 
species comprising the transverse Overhauser field. At 
low magnetic fields, this precession term leads to the col- 
lapse and revival effect. 



D. Stability of the revivals peaks against system 
inhomogeneities 

In this subsection we investigate how modifications of 
nuclear precession induced by the small spatial variations 
of the Knight fields and nuclear Zeeman couplings affect 
the electron spin echo signal, Eq. (|^^ . In particular, we 
focus on the stability of the revival peaks. 

Note that a necessary condition for the revivals to ap- 
pear is the clear separation of Larmor frequencies, such 
that ujki > SuJki,Aki for a{k) ^ a{l). For Sujki = Aki = 
0, i.e. for a homogeneous nuclear system with no intrin- 
sic nuclear spin dephasing, Eq. (|4T|) reduces to Eq. (f23| . 
where the right hand side exhibits revivals for values of 
the free evolution time r satisfying uJapr sa 47m, for 
any integer n. In this case the overall decay of the spin 
echo signal results solely from the spectral diffusion fac- 
tor e-^^/^^s")* in Eq. (|42l). 

Small but nonzero frequency differences Sujki and A^i 
can give rise to nuclear spin dephasing and cause addi- 
tional decay of the revival peak envelope as a function 
of T. We now analyze this decay and discuss its physical 
origin. Near the revival peaks LOkiT 4?™, and for early 
times r satisfying SuJkiT,AkiT <C 1, Eq. (|40|) simplifies to 



Tki\ peak 



ig* HBbjbkjujki + Su ki) (Amt)'^ - {SuJurf 
4 IB,, J 



[wki +5uJkiY - All' 
(43) 

For homo-nuclear terms with ujki ~ (i.e. for groups k 



rj. I _ ig* iiBhbkSi^ki 2 

4|i3oxt| 



a{k) = a(/), (44) 



while hetero-nuclear terms are given by 



Pkl I peak 



ig*liBbibk All - ^^li _2 
4|Bcxt| i^ki 



, a{k) ^ a{l). (45) 



Note that the hetero-nuclear terms and suppressed by 
the small ratios of or Hence, for 5ujkiT,AkiT < 
1, near the revival peaks the T- matrix is approximately 
block-diagonal with respect to the three species. Up to 
second order in UkiT, AkiT, we find that each block has a 
single pair of nonzero eigenvalues given by 



'^a Ipcak — ^ . , ^ , „ , T 



4iVrfLB, 



(46) 



where y^fJojJ) is the rms spread of Larmor angular fre- 
quencies of species a. Thus, we obtain a simple expres- 
sion for the echo envelope decay at the revival peaks: 



I peak 



n 



g* liBagnaA'i 



{5u 



Physically, Eq. 



(47) 



(|47|) describes decay of the revival enve- 
lope with timescale tJ^ = ( ^ an'!i\b'^'!\°' ) (^'^q)' which 
arises primarily from the intra-species spread of the Lar- 
mor frequencies, ^y{6uJ^. The effect of the off-diagonal 
matrix elements between different species is negligible. 
Additionally, the effect of the Knight field is also negligi- 
ble, due to the fact that the Knight field reverses its sign 
halfway through the evolution when the electron spin is 
fiipped by the 7r-pulse of the spin-echo sequence. 



E. Estimate of quantum corrections to the 
semi-classical results 



The semi-classical treatment is expected to be valid in 
the limit of a large number of participating nuclei, Nd- To 
better understand the validity of the approximation for 
finite Ndi in this section we provide a heuristic estimate 
for the deviation of the semi-classical expression for the 
single electron coherence function, Ccs, Eq. ([57)1 . from 
the quantum expression for C, given in Eq. psp . We 
define the quantum error as Ccs — C . 

We are interested in particular in the quantum error as- 
sociated with the semi-classical approximation to the dy- 
namics induced by the hyperfine and Zeeman couplings. 
Therefore we ignore the nuclear spin-diffusion contribu- 
tion to the decoherence, which is caused by the dipolar 
interaction. This is done by setting tsd = oo in Eqs. (|35p 
and p7p . We analyze the scaling of the quantum error as 
Nd is increased. However, while changing Nd, we wish to 
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keep fixed the rms values of the Overhauser field compo- 
nents, which determine the timescale for the decoherence 
of the electrons spins. Given Eq. (jG]), the rms value of 
the Overhauser field for each species scales as Aa/V^d, 
therefore we require Aa oc \/Nd. 

For large Nd, the leading contribution to the quantum 
error comes from the fact that for a given initial nuclear 
spin state — {^j, n^), the overlap of the final nu- 
clear spin states for different initial electron spin states, 
|*/+) = U±+ (t) I*) forcr = + and |*/_) = U±- (r) I*) 
for cr = —, is not unity, |( | ^P/- )| < 1. In other 
words, for every initial coherent state, the final state in- 
cludes quantum correlations between the electron and nu- 
clear spins, which are not captured by the semi-classical 
treatment. Thus, the electron-state-dependent modifica- 
tion to the nuclear evolution contributes an additional 
suppression of C which is not accounted for in Eq. (j36p . 

Up to leading order in the hyperfine coupling, the nu- 
clear spin evolution is dominated by Larmor precession 
due to the combination of the external magnetic field and 
the component of the Knight field parallel to the exter- 
nal field axis, z. Both contributions are included in Hq in 
Eq. ([2S1). However, these two effects alone will result in a 
perfect overlap of the final states, l^*/-!-) = l^'/-)- This is 
because the Larmor precession about the external field is 
independent of the electron spin state, while the preces- 
sion due to the Knight field is perfectly reversed halfway 
through the evolution due to the echo pulse which fiips 
the electron spin. Thus both effects were fully accounted 
for in the semi-classical treatment above. 

The quantum error results from the higher order terms 
in the hyperfine coupling, in particular from the trans- 
verse components of the Knight field, i.e. from the terms 
S'j"/,7 + S^I^ in the system's Hamiltonian, Eq. ([5|). 

These terms contribute to aHl{t) = ^^^^^^^j^^p^ in 
the reduced Hamiltonian, Eq. ([T^ . Although the trans- 
verse components of the Knight field are also reversed 
after the echo pulse, the fact that they do not commute 
with Hq means that difference in the final states 
and \^f-) may survive the echo. 

Below we first focus our discussion on a theoretical 
model which, in the semi-classical treatment, produces 
perfect revivals in the echo signal due to exact com- 
mensuration between Larmor periods of different nuclear 
species. We separate the discussion into two cases, in 
which the free evolution time r is either exactly on, or 
is away from, a revival peak. Then, at the end of the 
section we consider the quantum error at a revival peak 
which is not perfect, even within the semi-classical ap- 
proximation, due to a lack of commensuration between 
nuclear Larmor periods. 

For the initial state given above, which is a prod- 
uct of "giant spin" coherent states in each of the K groups 
of nuclear spins, we argue that, to leading order in Nd, 
the final states l^'/a) remain approximately tensor prod- 



ucts of coherent states, 

l*/.)«(8)|4,n,,/.) (48) 

k 

Within this picture, crH'^{t) causes the giant spin coher- 
ent states to rotate to new directions {nfej-o-}, which de- 
pend on the electron spin state cr = + or —. Other quan- 
tum effects such as coherent spin state squeezing due to 
the (/^)^ -I- {I^y^ terms in H'J^{t), which may stretch the 
coherent states anisotropically, or a build-up of quantum 
correlations between different giant spins, are also ex- 
pected. However, these effects enter only at higher order 
in I /Ik, because they are seeded by the small quantum 
fiuctuations of the spin components in the initial coherent 
state. 

Using the argument above, we find that away from the 
revival peaks the angle between the directions of the two 
final states of each giant spin fc, for a = + or — , scales 
as 1/Ik- Why is this so? First, according to Eq. ([5]), 
the Knight field acting on the nuclear spins, Ad.n^d, is 
smaller than the Overhauser field acting on the electron 
by a factor \/^/Nd ^ ^/Ik- Second, assuming that the 
timescale of the electron spin coherence collapse is com- 
parable to the timescale between revivals, the contribu- 
tion of H±_ {t) to the Overhauser field causes electron spin 
precession through an angle of order 27r. Over the same 
time interval, the corresponding part of the Knight field 
will cause the nuclear spins to rotate through an angle 
which is 1/Ik times smaller. 

The overlap | ( 5*/+ | VP/- ) | can thus be approximated 
to leading order by a product of overlaps between pairs 
of coherent states Y\k \{Ikii^kj+ \ Ik,i^k.f- )|, which are 
misaligned by an angle that scales down as 1 //fe . Further- 
more, because Ik is large, each coherent state \Ik,ii-kf±) 
is characterized by a Gaussian phase space distribution 
with angular width of order l/^/Tk- The reduction of the 
overlap, 1 — \{Ik,'n.kj+\Ik,'n.k.j-)\, therefore scales as 
l/A- ~ \J W2' therefore find that the overall reduc- 
tion of the total overlap, 1 — |( ^/+ | ^'/- )|, scales down 
at least as ~ ^/K^/Nd, as does the quantum error. 

However, exactly at a revival peak, where the condition 
of perfect nuclear precession commensurability is met, 
this leading contribution vanishes. Here, for a given gi- 
ant spin k, treating all other spins fc' 7^ fc as classical 
vectors, the precession and squeezing induced by II'[{t) 
is perfectly reversed after the 7r-pulse. Thus the quantum 
error at a revival peak is a higher order effect in 1/ Nd- 
Numerical simulations with just two giant spins indicate 
that the quantum error at the revival peaks may scale 
down even faster then l/Nd, but a more detailed analy- 
sis is a subject for further study. 

In more realistic cases, when the commensurability 
condition of all nuclear species cannot be exactly met, 
the quantum error near the quasi-revival peak scales like 
f{T)yjK^/Nj:, with the pre-factor /(t) becoming small 
as r approaches the approximate commensuration point. 
At the quasi-revival peak, /(r) is then dominated by the 
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classical effect of imperfect commensuration, captured by 
1 — Csc (t) . We estimate the error in the case of GaAs 
quantum dots with A^^ = (2 — 4) x 10^. For this esti- 
mate, we take the minimal number of groups of nuclear 
spins K which, within the semi-classical treatment, ac- 
curately produces the echo signal behavior near the re- 
vivals peaks, without changing significantly upon further 
refinement to more groups. For that purpose, note that 
only the spread of the Larmor frequencies for each nu- 
clear species, {5ujD, enters the expression for 1 — Csc (t) 
in Eq. (|Tf)). Consequently, it is sufficient to divide each 
of the three species into two collective groups. Group- 
in g the n uclei into K ^ 6 collective giant spins gives 
y^K^/Nd < 1/90. Therefore the quantum contribution 
to the imperfect revival is smaller than the semi-classical 
contribution by more than an order of magnitude. 

VII. DISCUSSION 

Equation (|42|) . together with Eq. (|40l). was used to fit 
the experimental echo signal data in Ref. O As stated 
in Sec. I VI Cl in general the matrix Tki and its eigenvalues 
{Am} are different for each dot due to variations in dot 
size and local environment. In particular, the distribu- 
tion of hyperfine couplings, {Afe}, and the number of nu- 
clei, Nd, depend on the distribution of electron density in 
the dot. Furthermore, the distribution of {SuJk} depends 
on, e.g. local electric field gradients which couple to the 
nuclear quadrupole moments. However, in practice we 
found that choosing the same parameters for Tki in the 
two dots produced a very good fit to the experimental 
data. Allowing, for example, different values of for 
the two dots did not significantly improve the fits. 

The spread of the Larmor frequencies y/ ((5cj^) that 
was found from fitting to the experimental resultai^ was 
equivalent to 3 Gauss effective spread of magnetic field, 
somewhat bigger than the values of NMR line-widths, 
typically about 1 Gauss equivalent spread, reported for 
bulk GaAs in the liter atur o'^^'^° . In addition to a random 
contribution to the local magnetic field coming from the 
dipolar interaction with neighboring nuclei, the nuclear 
spins in a GaAs quantum dot experience a quadrupolar 
splitting due to electric field gradients originating from 
the confined electrons, which we estimate to be at the 
order of a few Gauss (see Ref. [l^, supplementary mate- 
rials). These quadrupolar shifts may be responsible for 
the difference between the observed value and the NMR 
line-width. In the echo experiments which have been 
performed so far there is no way to distinguish between 
different origins of the apparent spread of nuclear Lar- 
mor frequencies (i.e. between decay due to nuclear spin 
fiip-flop, due to inhomogcneous broadening, or due to 
quadrupolar effects). 

Our model does not directly include the dipole-dipole- 
induced temporal decay of local nuclear spin correlations 



(/,\(0)/j((i)). However, such decay is accounted for phe- 
nomenologically by a time-independent spread of site- 
dependent Larmor precession frequencies. Because the 
current Hahn echo experiments can't distinguish between 
these intrinsic and extrinsic nuclear spin dcphasing pro- 
cesses, our approximation of accounting for these effects 
by a random static, disordered, Zceman field is reason- 
able. According to Eq. (|17| . the affect of this random 
field is to cause an additional decay to the echo signal 
with a timescale t±. This decay becomes dominant at 
low magnetic fields, as t± becomes shorter then the de- 
cay time associated with the spectral diffusion, tsd- 

In addition, note that the spectral diffusion decay time 
TSD was found experimentally to be independent of the 
magnetic field strength. Other experiments^ have also 
suggested that the external magnetic field does not signif- 
icantly influence spin-diffusion at the relevant field range, 
above 20 mT. This provides further justification for ne- 
glecting the commutators [H± , Hd] in section IVIl which, 
a priori, could introduce such a dependence. 

In conclusion, we have provided physical justification 
for treating the Overhauser fields in a GaAs double quan- 
tum dot system as classical vector variables, based on 
the relative smallness of the commutation relations due 
to the large number of spins in each dot. For the sim- 
ple model in section |Vl in which we ignored the nuclear 
dipole-dipole interaction and the Knight shift of the nu- 
clear Larmor frequencies, we demonstrated the equiva- 
lence of the semi-classical treatment and the quantum 
diagrammatic summation of Ref. [32l . The semi-classical 
treatment quantitatively captures the observed phenom- 
ena of monotonic decay of the echo signal in strong mag- 
netic fields, and the collapses and revivals of the echo 
signal in weaker magnetic fields. The overall effect of the 
Knight field on the spin-echo is found to be negligible. 
This fact is parameter dependent. However note that the 
Knight field is reversed by the echo pulse. As a result. If 
one treats the Knight field contribution to the phase in 
Eqs. and P5|) as a perturbation, it vanishes in the 
leading order in r, near the revivals peak. This is unlike 
the contribution of competing process, i.e. the spread of 
the Larmor frequencies. This means that, for the typical 
parameters in GaAs quantum dots, after averaging over 
all nuclear spin states, both the collapse-revivals effect 
and the overall envelope decay can be understood simply 
as arising from averaging over the initial conditions of 
the Larmor precession of the nuclear spins. 
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